Plant Molecular Biology

In wheat, weak seed dormancy (SD) is related to an increased tendency for pre-harvest sprouting (PHS), which reduces yield and quality. However, the molecular mechanism underlying SD remains elusive. Here, we identified a wheat R2R3-MYB transcription factor (TaMYB83-7B) related to SD. Expression analysis showed that TaMYB83-7B was highly expressed in wheat seeds, and was more highly expressed in strong-dormancy varieties than in weak-dormancy varieties. Sequence and association analysis indicated that T/C mutations at -907 bp and -1133 bp in the TaMYB83-7B promoter were significantly associated with wheat SD, with C at both sites related to strong dormancy. Dual-luciferase reporter assays demonstrated that the transcriptional activity of the TaMYB83-7B promoter was significantly higher in strong-dormancy varieties than in weak-dormancy varieties. Further analyses indicated that TaMYB83-7B functions as a transcriptional inhibitor. Germination experiments revealed that overexpression of TaMYB83-7B significantly enhanced SD, while its loss-of-function reduced SD. Finally, TaMYB83-7B was found to regulate SD by influencing the balance between abscisic acid (ABA) and gibberellin (GA) in wheat seeds. Overall, the results of this study enhance our understanding of the complex regulatory mechanism underlying SD, and provide gene targets and molecular markers for the genetic improvement of PHS resistance in wheat.

RAP2.6, an AP2/ERF transcription factor (TF), regulates plant stress responses; however, its role in floral transition remains unexplored. Here, we evaluated RAP2.6s role in flowering and the associated transcriptional changes in Arabidopsis thaliana under long-day conditions. RAP2.6-overexpressing line showed early flowering with fewer rosette leaves, whereas rap2.6-1 mutant flowered later, had more rosette leaves, and higher expression of the floral repressor FLOWERING LOCUS C (FLC). Early flowering in the overexpressing line was accompanied by transcriptional activation of the floral integrators GIGANTEA (GI), FLOWERING LOCUS T (FT), and COSTANS (CO), potentially through RAP2.6 interaction with GCC/DRE cis-regulatory elements. RAP2.6-mediated floral transition depended on nitric oxide (NO), with flowering time largely varying based on NO bioactivity. RAP2.6 was found to be a downstream regulator of Arabidopsis S-NITROSOGLUTATHIONE REDUCTASE 1 (GSNOR1) in controlling S-nitrosothiol (SNO) levels, flowering time, and silique formation. The NITRIC OXIDE-ASSOCIATED 1 (NOA1)-dependent reduction in NO levels abolished early flowering in 35S::RAP2.6 plants without affecting silique formation. Furthermore, enhanced cytokinin sensitivity and upregulation of cytokinin biosynthetic genes suggest cytokinin involvement in RAP2.6-mediated flowering. Together, these findings highlight the crucial role of RAP2.6 in regulating flowering time by integrating redox and hormonal signaling to coordinate reproductive development in A. thaliana.

Obtaining tomato plants with firm and intact fruit is one of the main goals in tomato breeding programs. Achieving these goals through conventional breeding is time-consuming and can lead to the loss of unwanted traits. In other hand, consumers are concerned about the presence of transgenic elements in plants acquired through RNA interference. The use of CRISPR/Cas9 technology has made it possible to overcome the above-mentioned shortcomings. In this study, the {beta}-D-N-acetylhexosaminidase ({beta}-hex) gene, which is involved in tomato fruit ripening, was knocked out using CRISPR/Cas9. In the resulting mutant plant genome, an indel mutation was found in exons 1 and 2 of the {beta}-hex gene. Plants with a mutation in their genome were observed to have increased fruit firmness and shelf life compared to control plants without affecting fruit quality.

We recently reported the identification of endogenous viral elements (EVEs) originating from the Caulimoviridae family within the alfalfa (Medicago sativa L.) genome. Our subsequent identification of ubiquitous rhabdoviral elements in infected and healthy alfalfa tissues by high throughput sequencing prompted us to suggest that the alfalfa genome might be populated with integrated rhabdoviruses as well. Bioinformatics analysis using 26 publicly available alfalfa genomes proved the suggestion accurate. We found multiple non-retroviral segments of the Rhabdoviridae family belonging to the genera Betanucleorhabdovirus and Betacytorhabdovirus that appeared to be stable constituents of the host genome. In that capacity they could potentially acquire functional roles in alfalfas development and response to environmental stresses. We believe this study reveals the first documented case of rhabdoviruses integrated into the alfalfa genome.

Salinity is a major crop production constraint in dry pea (Pisum sativum L.), making the development of salt-tolerant varieties essential to improve crop productivity and land-use efficiency. The genetic mechanisms of salt tolerance in dry pea is largely unknown, and research on salt-tolerant genes is limited. In this study, we established comprehensive genomic and transcriptomic resources, along with a robust screening protocol, to dissect the genetic basis of salinity tolerance using two germplasm sets: the USDA pea diversity panel, consisting of approximately 200 globally sourced accessions, and a set of 300 modern elite lines from the NDSU Pulse Crops Breeding Program. Genetic variation for the salinity response was assessed based on ten phenotypic traits, with root dry weight, shoot dry weight, and specific root length identified as key indicators based on their heritability. Genome-wide association mapping uncovered significant genomic regions and several candidate genes linked to salt stress, with the strongest association found on chromosome 6. Overlapping QTL signals across traits suggest a shared genetic architecture underlying salinity tolerance. Field-based transcriptomic analysis further identified five putative genes involved in salinity response conserved across multiple crop species. Notably, Psat5g000800, encoding a glycosyl hydrolase gene, was markedly upregulated under salinity stress. These findings highlight the complex, multi-gene regulatory nature of salinity tolerance in dry pea and underscore the importance of functional validation of candidate genes. This study provides key insights and practical tools to support breeding efforts aimed at improving salt tolerance in dry pea.

Tomato wild relatives are valuable genetic resources for trait discovery and understanding the genetic basis of fruit metabolism and quality. Yet, only a fraction of naturally occurring variation has been exploited. Here, we performed metabolite profiling of two large Backcross Inbred Line populations derived from crosses between the wild species S. pennellii accession LA5240 (Lost) and cultivated genotypes LEA (determinate) and TOP (indeterminate), including [~]1400 and [~]500 lines, respectively. High-resolution mapping identified enormous metabolic quantitative trait loci (mQTL), including a new locus on chromosome 12 associated with fruit sucrose accumulation that harbours INVERTASE INHIBITOR 3 (SlINVINH3) protein. Comparative analysis indicated that SlINVINH3 is highly expressed in wild S. pennellii 0716 fruit, whereas a six-amino acid deletion is present in its coding sequence compared with S.pennellii LA5240 and S. lycopersicum. We further demonstrated that in SlINVINH3-overexpressing tomato plants, only the S. pennellii LA5240 allele led to increased sucrose, accompanied by reduced fructose and glucose levels. Furthermore, the large population size enabled us to assess the epistatic interactions, with approximately 40% of interactions being more-than-additive and 60% less-than-additive. Our results demonstrate the power of permanent exotic populations to reveal hidden metabolic diversity and provide an approach for improving fruit quality through targeted breeding and metabolic engineering.

Magnaporthe oryzae, the rice blast fungus, plays a role as a model organism for molecular plant-microbe interaction research. Studies on the pathogenic mechanism of this fungus revealed many genes involved in signaling pathways. As multi-omics data are being available, genomic-level researches have been conducted to uncover the underlying biological processes during the pathogenesis of M. oryzae. Identifying the genome-wide protein-protein interaction (PPI) network is one of the omics-level approaches, which helps to understand signaling and regulatory pathways. However, existing biological network resources of M. oryzae are not sufficient to decipher pathogenesis mechanisms due to the abundance of false positives/negatives. In this study, a reliable PPI network database of M. oryzae, MagNet, was constructed with three methods, including homology-based Interolog search, co-expression network construction, and domain-domain interaction (DDI)-based prediction. With three approaches altogether, the pan-network with 5,600,976 interactions was generated, including 217,531 highly confident interactions supported by all three methods. Experimental data on M. oryzae PPIs supported that our PPI network can predict PPIs with higher accuracy compared to the previously constructed databases. MagNet would provide integrated biological network data, which can help to understand the molecular mechanisms of the rice blast fungus. The PPI data can be accessed via https:/magnet.scnu.ac.kr.

Here, we present a comprehensive multiomics analysis of anthocyanin biosynthesis in Rubus armeniacus, known for its dark fruits. A phased genome sequence of the tetraploid blackberry was generated, achieving an N50 of 34 Mb with an assembly size of 1.2 Gbp based on Oxford Nanopore Technology sequencing (ONT). The BUSCO score for the total assembly shows a high completeness of 99.1%. The assembly was separated into 4 pseudohaplophases, with the pseudohaplophase A representing the R. armeniacus genome in 7 chromosome scale contigs, with an N50 of 46 Mbp and 98.8% conserved BUSCO genes. A total of 118,183 protein coding genes were annotated within the genome assembly and all relevant genes encoding enzymes and transcriptional regulators of the anthocyanin biosynthesis pathway were identified within each pseudohaplophase. To further understand the underlying cause of dark pigmentation, the gene expression was analysed during different stages of berry development revealing a strong induction of anthocyanin biosynthesis genes including the anthocyanin activating subgroup 6 MYB transcriptions during the berry ripening process. Further, a quantification of cyanidin-3-O-glucoside in methanolic berry extract, utilizing a UHPLC-HRAM-MS analysis, revealed an approximately 500-fold increase of cyanidin-3-O-glucoside from green to black fruit, indicating that dark pigmentation in R. armeniacus results from high anthocyanin accumulation. Significance statementThis study provides a multiomics analysis of the dark pigmentation of Rubus armeniacus, including a high quality phased assembly and an in-depth analysis of the anthocyanin biosynthesis pathway. A transcriptional and metabolomic analysis revealed that dark berry pigmentation is caused by a high accumulation of cyanidin-3-O-glucoside during fruit ripening.

Wild bees face declines, and forests may serve as critical habitats for pollinators. However, how forest composition and the associated floral environment shape pollen provisioning and resource partitioning among cavity-nesting bees remains poorly understood. Here, we leveraged BEF-China, a large-scale subtropical forest biodiversity experiment with experimentally controlled plant (tree and shrub) communities, to investigate how forest composition and spatial context shape pollen provisioning, resource partitioning, and reproductive success of cavity-nesting bees. We used DNA metabarcoding to analyze floral composition of pollen provisioned by five cavity-nesting bee species, with samples collected from BEF-China across three years (2022- 2024). By comparing pollen taxonomic composition from whole-nest pooled samples and individual brood-cell samples with the experimentally planted species pool, we characterized dietary patterns and temporal dynamics of five bee species. Bees primarily relied on floral resources from the surrounding landscape, with planted trees providing essential but temporally restricted pollen supplements during specific phenological stages. Co-occurring bee species exhibited staggered nesting phenology and distinct dietary preferences for different plant families, with fine-scale resource differentiation even during periods of phenological overlap. Our results suggest that managed forests support cavity-nesting bees by providing critical woody floral resources during specific phenological gaps and offering stable nesting environments. To mitigate pollinator declines, forest management should prioritize maintaining diverse, phenologically complementary flowering vegetation within and surrounding forest stands. This ensures temporal continuity of pollen availability throughout the nesting season, which is particularly crucial for restoring pollinator services in simplified forest landscapes.

RNA-seq datasets from medicinal yews are crucial for studying paclitaxel biosynthesis. However, cross-study data analyses are hindered by pronounced batch effects. Here, we compiled 45 RNA-seq samples from three studies across four tissues (bark, leaf, root, stem) and assessed 35 preprocessing pipelines combining six normalization strategies with five batch-effect correction approaches. Unsupervised clustering (HCA, k-means, Grade-of-Membership), evaluated using Jaccard and Adjusted Rand indices, revealed significant variability in batch effect removal. Supervised classification of tissue and project labels (Random Forest and linear/radial SVM) demonstrated improved accuracy in tissue type prediction, highlighting the effectiveness of correction methods. The processed data facilitated the identification of 189 putative ABC transporters across samples, six of which showing a strong correlation to the gene encoding 10-deacetylbaccatin-III-10{beta}-O-acetyltransferase, a key biosynthetic enzyme in the taxol pathway. High expression levels in leaf and bark further support their role in taxane intermediates trafficking in taxol biosynthesis. Structural analysis and molecular docking further supported the selection of these candidates, and the agreement between transcriptomic ranking and docking-based prioritization suggests that these transporters may participate in taxane intermediate recognition, trafficking, or export. These findings demonstrate the importance of normalization and batch effect correction in RNA-seq analysis to advance gene discovery in Taxus species and, more broadly, in plant research. Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=152 SRC="FIGDIR/small/723993v1_ufig1.gif" ALT="Figure 1"> View larger version (54K): org.highwire.dtl.DTLVardef@1469162org.highwire.dtl.DTLVardef@1f2c4deorg.highwire.dtl.DTLVardef@15ad821org.highwire.dtl.DTLVardef@123676d_HPS_FORMAT_FIGEXP M_FIG C_FIG

The production of reactive oxygen species (ROS) in response to cadmium (Cd) has been extensively studied, demonstrating that they play a key role in the plants response to this heavy metal. While the role of enzymes like RBOHs has been thoroughly studied, the function of other ROS-producing enzymes, such as peroxisomal glycolate oxidase (GOX), remains largely overlooked. Peroxisomal GOX is a core metabolic enzyme of the photorespiratory pathway occurring in chloroplasts, mitochondria and peroxisomes. Using Arabidopsis (Arabidopsis thaliana) mutants lacking the main peroxisomal GOX genes, GOX1 (gox1-1) and GOX2 (gox2-1) we explored their function in plant response to Cd. Although photosynthetic capacity appears to be affected to the same extent in both mutants under control and Cd stress conditions, GOX2 seems to play a greater role in ROS production in response to the metal. Transcriptomic analyses on WT and gox2-1 pointed to the mitochondrial electron transport chain (mETC) as a target of Cd stress. We further investigated the individual GOX1 and GOX2 functions in mETC regulation and redox state. Although oxidative ratio of mitochondria was higher in both mutants, it was more pronounced in the absence of GOX1. Furthermore, the mETC is affected in both mutants but the regulation of its components differs in each mutant. These results point out the different functions of the two photorespiratory GOX isoforms in Arabidopsis, leading to a better understanding of the photorespiratory pathway.

In sunflower (Helianthus annuus L.), the composition of fatty acids in the seeds, primarily oleic, linoleic, stearic and palmitic acid, is of utmost importance for oil quality. Despite this, the genetic basis of this trait and its interaction with the environment is poorly understood. Understanding this interaction is critical to improvement of sunflower within the context of climate change. In this work, we incorporated fatty acid composition measurements from the sunflower SAM population and eight environments across an extensive geographic cline into GWAS. The SAM panel consists of 287 varieties representing approximately 90% of sunflower diversity, for which 2.2 million high-quality SNPs with a MAF > 5% are available. For increased power, multivariate GWAS was performed with four different inputs: (i) mean fatty acid composition within each environment, (ii) mean fatty acid composition within each environment omitting high oleic varieties, (iii) trait stability within environments quantified by standard errors among replicate samples ( stability) and (iv) Eberhart and Russells {beta} which quantifies trait stabilities across environments ({beta} stability). All four analyses yielded highly significantly associated SNPs. We found that high oleic varieties exhibited high {beta} trait stability, resulting in substantial overlap in markers between analyses (i) and (iv), with signals being fairly consistent between environments in analysis (i). For analyses (ii) and (iii), significant markers tended to vary between trials. For significant SNPs across all analyses, 147 candidate genes were identified, including promising candidates such as 15 fatty acid metabolism genes, 6 heat shock proteins and 22 transcription factors. Lastly, a large introgression consisting of two flanking inverted sequences on Chromosome 5 was found to coincide with stability in the Georgia trial, suggesting a role in FA composition stability under high heat conditions.

Chilling stress induces photosystem I (PSI) photoinhibition in chilling-sensitive cucumber, in which insufficient activity of the chloroplast NADH dehydrogenase-like complex (NDH) leads to PSI over-reduction and damage. However, it is not yet clear whether these findings can be generalized to other species or what the molecular mechanism underlying impaired NDH function is. In this study, we first examined whether NDH is essential for PSI protection under chilling stress using an NDH-deficient rice mutant. Compared with wild-type plants, the NDH-deficient mutant exhibited enhanced PSI over-reduction and pronounced PSI photoinhibition under chilling stress. In contrast, rice plants expressing flavodiiron protein (FLV), which functions as an alternative electron acceptor downstream of PSI, did not exhibit PSI photoinhibition under chilling stress, demonstrating that electron sink capacity of NDH is important for PSI protection under chilling stress. Furthermore, analysis of the factors responsible for NDH dysfunction under chilling stress in cucumber revealed that chilling stress destabilizes the PSI-NDH supercomplex, leading to NDH monomerization and a consequent loss of NDH activity. This NDH monomerization is likely attributable to chilling-induced damage to the light-harvesting complex Lhca, which mediates the association between PSI and NDH. Together, these results indicate that NDH is essential for protecting PSI from photoinhibition under chilling stress in both rice and cucumber, and that chilling-induced destabilization of the PSI-NDH supercomplex represents a key molecular mechanism underlying PSI over-reduction and photoinhibition.

Chlamydomonas reinhardtii is a model green alga extensively used to study photosynthesis and cilia using molecular biology and genetics. Electroporation is a very common technique to transform DNA into the nuclear genome, which is essential to generate mutant collections and express transgenes. Here, we describe a simple, fast, and efficient protocol to transform strains with an intact cell wall. It achieves a good transformation efficiency without cell wall digestion or use of commercial kits and is compatible with the widely available Gene Pulser electroporation system. Key featuresO_LIHigh transformation efficiency of Chlamydomonas reinhardtii strains with an intact cell wall. C_LIO_LIFaster than currently available electroporation protocols. C_LI

Hemp (Cannabis sativa) produces a wide array of medicinally significant compounds, including cannabidiol (CBD). These compounds are predominantly synthesized in female hemp inflorescences. The proposed research utilizes next-generation sequencing-based transcriptome analysis using a 3{square}-end-directed approach to identify differentially expressed genes between male and female hemp plants at the early vegetative stage. 886 differentially expressed genes (DEGs) were identified, a majority of which were upregulated in males compared to females. We hypothesized that alternative RNA processing contributes to sex-specific gene expression. To this end, 932 genes were identified that exhibited significant changes in poly(A) site usage when comparing males and females. These genes were much more likely to be differentially expressed, supportive of this hypothesis. Males tend to have longer 3 UTRs with canonical motifs found in the Near-Upstream Elements (NUE), compared to the shorter 3 UTRs in females, which have A-rich motifs near the cleavage site. This suggests that polyadenylation remodels hemp mRNAs with distal poly(A) sites being preferred in males. To further investigate when this sex-specific gene expression program is established, RNA was isolated from plants at various developmental stages, such as developing seeds, four-day-old seedlings, and different developmental stages up to four weeks after sowing. Diagnostic male-specific genes were analyzed using RT/PCR. The results indicate that sex-specific gene expression is not evident in seeds but rather is set during or after germination. SignificanceO_LIHemp males tend to have longer 3 UTRs with canonical motifs found in the Near-Upstream Elements (NUE), compared to the shorter 3 UTRs in females, which have A-rich motifs near the cleavage site. C_LIO_LIThe sex-specific gene expression program is not yet established in mature seed but is set in the time between germination and 4 days of growth. C_LI

Pennycress (Thlaspi arvense L.) is an intermediate winter oilseed crop that has only recently been domesticated for agronomic use. Improving agronomic traits requires sources of genetic variation, and mutagenesis is frequently used to help overcome the limitations of natural populations. We investigate the impact of Ethyl methanesulfonate (EMS) on genetically effective cells (GECs) to characterize the intra-individual genetic variation of EMS mutagenesis in pennycress. We identified that pennycress contains at least 4 GECs which, when treated with EMS, create unique mutations across different branches within the same individual plant. We then propagated the M2 plants for whole genome sequencing, providing extensive characterization of the EMS mutation profile and developing a gene index as a resource for future reverse genetic screenings. Article SummaryPennycress is an emerging winter oil seed crop in the American Midwest. Domestication efforts have advanced rapidly through a combination of genetic techniques. One of the most successful methods has been the use of a mutant gene index, a large collection of pennycress seed where new genetic variation has been created through Ethyl methanesulfonate (EMS). EMS mutations are not uniform however, and a single treated seed can have wide genetic variation within the resulting plant. We investigate the role of genetically effective cells on EMS variation, and present the full EMS population as a resource for further pennycress domestication efforts.

Plant receptor-like kinases (RLKs) are involved in diverse processes, ranging from growth and reproduction to interactions with microbes. Variation in the extracellular domains delineates several RLKs subfamilies, including the malectin-like domain leucine-rich repeat receptor-like kinases (MLD-LRR-RLKs). Symbiosis Receptor-like Kinase (SymRK) is the prototypical member of MLD-LRR-RLKs and is required for microbial accommodation in host roots during root endosymbiosis. Yet, comparative phylogenetic analysis of SymRK orthologs in the broader context of MLD-LRR-RLK subfamily evolution remains limited. In this study, we examined the inventory, phylogeny and clade-specific evolutionary and transcriptional characteristics of this receptor group. SymRK and its closest homologs are present in most land plant lineages and group into four major clades and six additional species-specific clades. These clades can be distinguished by their evolutionary characteristics as either conserved with reduced gene copy number changes (including SymRK) or expanded and diversified, as observed in clade IV. Clade IV dynamics are largely driven by tandem gene duplications, which often arise within gene clusters. We further analysed the evolutionary characteristics of MLD-LRR-RLKs at the population level in Arabidopsis thaliana accessions. We found that some genes are conserved across accessions and are therefore likely to be functionally important, whereas a subset of genes, often located within tandem clusters, are highly diverse and likely contribute to accession-specific adaptations. Finally, most MLD-LRR-RLKs in the A. thaliana Col-0 accession are expressed in roots and respond broadly to biotic stimuli at the transcriptional level. Notably, clustered genes frequently exhibited divergent expression profiles, suggesting transcriptional diversification. Together, we revealed two contrasting evolutionary characteristics among members of the MLD-LRR-RLK subfamily, potentially associated with their functions in plants.

Prenylated isoflavonoids are widely distributed specialized metabolites within the Fabaceae and contribute to various characteristic biological activities for both plants and humans. Several aromatic prenyltransferases (PTs) have been identified in Glycyrrhiza species, which are the most widely consumed crude drugs in traditional Chinese medicine. However, these enzymes do not sufficiently explain the structural diversity of prenylated flavonoids produced in the Glycyrrhiza genus. To identify additional novel PTs, we used elicited cultured Glycyrrhiza glabra roots as source material, in which elicitor treatment of cultured roots increased the accumulation of multiple prenylated flavonoids. To identify the responsible enzyme, PT candidates were screened using G. uralensis transcriptomes, currently the sole publicly available transcriptomic resource within the genus, and a homolog designated GgBSPT1 (BSPT; a broad-substrate prenyltransferase) was subsequently isolated from elicited cultured G. glabra roots. GgBSPT1 differed from previously identified Glycyrrhiza PTs in both amino acid sequence and enzymatic properties. GgBSPT1 catalyzed 3'-prenylation of isoliquiritigenin and 6-prenylation of five flavonoids, i.e., this PT displayed broad substrate acceptance across 20 distinct flavonoid structures. Overall, elicited cultured G. glabra roots enabled the identification of a previously unrecognized PT that is functionally distinct from earlier reported Glycyrrhiza PTs. This study provides a new insight into the metabolic plasticity of Glycyrrhiza species and expands the enzymatic toolkit for future metabolic engineering of prenylated phytochemicals by the unusually broad substrate specificity of GgBSPT1. Main conclusionUsing cultured Glycyrrhiza glabra roots, we identified a new prenyltransferase involved in the formation of a variety of flavonoids, thereby revealing novel prenylated isoflavonoid pathways in licorice.

Pangenome-level gene family identification often applies sequence similarity clustering without phylogenetic or synteny information, which risks biologically misleading evolutionary inferences. Using five transcription factor families (bHLH, MYB, NAC, WRKY, MADS-box) across 401 rice pangenome accessions, we compared clustering strategies: OrthoFinder alone, cd-hit alone, MMseqs2 alone, and OrthoFinder-informed refinement by cd-hit or MMseqs2. Methods solely based on sequence similarity merged distinct orthogroups and generated fewer orthogroups than approaches incorporating graph-based orthology. Conflicting cluster assignments, measured against OrthoFinder, varied strongly among families, from approximately 14% in MADS-box to approximately 57% in MYB, and were associated with protein length differences. Core, shell, and cloud gene classifications shifted substantially depending on the method, especially in MYB, NAC, and WRKY families. Critically, Ka/Ks distributions for core genes were highly method-sensitive, with orthology-aware methods yielding more convergent and less variable estimates of selective pressure, whereas noncore gene estimates remained robust. These findings demonstrate that neglecting graph-based orthogroup inference inflates methodological artifacts. We recommend a two-step strategy: initial graph-based orthogroup delineation followed by sequence similarity refinement to balance evolutionary accuracy and resolution in pangenome-scale gene family studies.

Indica rice, being a tropical crop, is highly sensitive to cold temperature. Cold stress affects vegetative growth, photosynthetic efficiency, along with reproductive features. Genetic resource screening in diverse landraces is an approach for identifying cold-tolerant traits. Here, we have characterised a boro germplasm, CB1, with an efficient germination rate and growth vigour when treated at chilling temperatures. CB1 seedlings show a higher survival rate compared to IR36 when subjected to prolonged chilling stress. Biochemical analyses indicated efficient ROS modulation, higher chlorophyll content, enhanced photosystem II efficiency and unique stomatal traits, leading to higher relative water content in CB1 plants during stress and recovery. Transcriptome analysis supported upregulation of chlorophyll biosynthesis, photosystem, & light harvesting complex and ROS scavenger genes in CB1 seedlings. Interestingly, high D1 protein turnover in CB1 promotes damage-repair of PSII for efficient photosynthesis. Furthermore, key transcription factors for stomatal development and expression of photosynthetic genes were upregulated in CB1 during stress recovery. Notably, higher expression of OsGLK1 and enrichment of GLK1 targets were observed in CB1 plants during chilling stress and recovery. Taken together, our results suggested that CB1 plants exhibit cold tolerance by modulating photosynthesis efficiency and stomatal behavior for better adaptability and survival against chilling temperature. HIGHLIGHTSThe efficient photosynthetic recovery, active ROS scavenging system and maintenance of water content through regulating stomatal traits, enhance the survival of indica germplasm CB1 against chilling stress.